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ABSTRACT 

The recent discovery of dust-correlated diffuse microwave emission has prompted two rival explana- 
tions: free-free emission and spinning dust grains. We present new detections of this component at 10 
and 15 GHz by the switched-beam Tenerife experiment. The data show a turnover in the spectrum and 
thereby supports the spinning dust hypothesis. We also present a significant detection of synchrotron ra- 
diation at 10 GHz, useful for normalizing foreground contamination of CMB experiments at high-galactic 
latitudes. 

Subject headings: cosmic microwave background - diffuse radiation - radiation mechanisms: thermal 
and non-thermal - methods: data analysis 



1. INTRODUCTION 

Understanding the diffuse microwave emission from the 
Galaxy is crucial for doing cosmology with Cosmic Mi- 
crowave Background (CMB) anisotropics. Although three 
components of Galactic emission have been firmly iden- 
tified (synchrotron and free-free radiation, and thermal 
emission from dust particles), it is important to better 
quantify their frequency dependence and spatial distribu- 
tion (see, e.g., Kogut et al. 1996a; Tegmark et al. 1999 and 
references therein). 

Cross-correlations of CMB data with far-IR maps have 
shown the existence of a microwave emission component 
whose spatial distribution is traced by these maps (Kogut 
et al. 1996b; Lim et al. 1996; de Oliveira-Costa et al. 1997; 
Leitch et al. 1997; de Oliveira-Costa et al. 1998, hereafter 
dOC98). Although this emission component has a spec- 
tral index suggestive of free- free emission (Kogut 1999), 
the correlations between Ha and dust maps have been 
found to be marginal (McCullough 1997; Kogut 1997). 
The source of this correlated emission is therefore an open 
question. Recent work suggests that it originates from 
spinning dust grain emission (Draine and Lazarian 1998), 
which should have a spectral index of —3.3 < {3 sp i n < — 4 
between 19 and 53 GHz. Since spinning dust grains have a 
predicted turndown in their emission at frequencies below 
20 GHz, a cross-correlation analysis with lower frequency 
measurements may help to discriminate between free-free 
and spinning dust emission models. 

The purpose of this letter is to evaluate the Galactic 
contribution in the Tenerife 10 and 15 GHz data sets by 
cross-correlating them with the DIRBE dust maps and 
with the Haslam and Reich and Reich synchrotron maps. 
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Our analysis is based on the Tenerife measurements made 
up until the end of 1997 (Gutierrez et al. 1999, hereafter 
G99). The Tenerife switched-beam experiment is a sky 
survey between 0° < RA < 360° and 30° < 6 < 45° 
(Figure 1 shows a sample strip at S = 32°) carried out at 
an angular resolution of 5.1° FWHM with an instrument 
that uses a double-differencing technique. Data was taken 
at frequencies 10 and 15 GHz (we omit the 33 GHz data 
here due to lack of low latitude sky coverage), which are 
treated separately in order to gain additional frequency 
information on Galactic emission. 
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Fig. 1 — Scan of the Tenerife region 0° < RA < 360° for 
8 = 32°. The thin curve shows Tenerife 15 GHz data, while 
the thicker curve represents DIRBE 100/im convolved with 
the Tenerife triple-beam, rescaled for best fit. The spikes at 
RA~ 80° and 300° correspond to the Galactic plane crossings. 



2. METHOD 

The Tenerife data consists of N = 2880 pixels at each 
frequency with double-differenced sky temperatures and 
noise rii. We assume that this data is a linear superposition 
of CMB fluctuations x l CMB and M Galactic components 
whose angular distributions are traced in part by exter- 
nal foreground templates. Writing these contributions as 
A-dimensional vectors, we obtain 



y = Xa 4- x C mb 



(1) 
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where X is an TV x M matrix whose rows contain the 
various foreground templates convolved with the Tenerife 
triple-beam (i.e., Xy would be the i th observation if the 
sky had looked like the j th foreground template), and a 
is a vector of size M that gives the levels at which these 
foreground templates are present in the Tenerife data. 

We treat n and xcmb as uncorrelated random vectors 
with zero mean, and the matrix X as constant; thus the 
data covariance matrix is given by 



C = (yy T ) - (y)(y T ) = (x CMB x£ MS ) + (nn 5 



(2) 



The Tenerife noise covariance matrix (nn T ) is, to good 
approximation, diagonal (see Gutierrez de la Cruz et al. 
1995), while the CMB covariance matrix has the form 
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where the vector w = (—0.5,1,-0.5) gives the weights 
of the three lobes in the Tenerife triple beam, f jo is the 
direction that the central lobe pointed to during the i th 
observation, and r^ii is at the same declination but off- 
set by ±8.1° on the sky. We take the CMB correlation 
function to be 



c(r ik 



c(cos 9) 
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B$C e Pi{cos9) (4) 



for a flat power spectrum Ce cx l/C(£ + 1) normalized to 
a Q flat = (5C2/471-) 1 / 2 = 20/iK (G99). Finally, we model 
the Tenerife beam as a Fisher function (Fisher et al. 1987) 
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with FWHM= V8 In 2a = 5.1°, which gives a window 
function Bi w e~" i ( t + 1 )/' 2 . The Fisher function is locally 
Gaussian and integrates to unity over the unit sphere. 

Since our goal is to measure a, both n and xcmb act 
as unwanted noise in equation ([!]). Minimizing \ 2 = 
(y — Xa) T C _1 (y — Xa) yields the minimum- variance es- 
timate of a, 
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a = [x T cr 1 X] 

with covariance matrix 

S=(a 2 )-(a) 2 =[X r C- 1 X]- 1 . (7) 



The error bars on individual correlations are there- 
fore Acii = . This includes the effect of chance 

alignments between the CMB and the various template 
maps, since the CMB anisotropy term is incorporated in 
(xcmb^cmb)- 

3. DATA ANALYSIS & RESULTS 

Due to the double-differencing technique, the Tenerife 
data are insensitive to the monopole (£ = 0) and dipole 
(£ — 1). Accordingly, when we convolve the template 
maps with the Tenerife triple beam function, we are re- 
moving the mean of the templates as well as large angular 



scale structure. As a consequence, our results depend pre- 
dominantly on the small scale intensity variations in the 
templates (i ~ 10 — 30) and are insensitive to the zero 
levels and gradients in the CMB and the template maps. 



Table 1. - Correlations for 10 and 15 GHz data. 
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The DIRBE and Haslam correlations listed in this table cor- 
respond to joint 100/ini— Has fits, whereas the R&R numbers 
correspond to a joint 100/iin— R&R fit. 

(b) a has units /xK (MJy/sr) -1 for the 100/im template, ^iK/K 
for the Has template and fiK/mK for the R&R template. We 
use antenna temperature throughout. 
< c > AT = (a ± Aa)crG a ;. 



3.1. Correlations & Variances 

We cross-correlate the Tenerife data with two differ- 
ent synchrotron templates: the 408 MHz survey (Haslam 
et al. 1981) and the 1420 MHz survey (Reich and Reich 
1988), hereafter Has and R&R, respectively. To study dust 
and/or free- free emission, we cross-correlate the Tenerife 
data with three Diffuse Infrared Background Experiment 
(DIRBE) sky maps at wavelengths 100, 140 and 240 ^m 
(Boggess et al. 1992). The extent of point source contam- 
ination in the Tenerife data is discussed and estimated in 
G99, and therefore will not be addressed in this Letter. 
In practice, we just remove the estimated point sources 
contribution before calculating the correlations. For def- 
initeness, we use the DIRBE 100/im channel when plac- 
ing all limits below since it is the least noisy of the three 
DIRBE channels, and the Haslam map since it is the syn- 
chrotron template at lowest frequency. 

Table 1 shows the coefficients a and the corresponding 
fluctuations in antenna temperature in the Tenerife data 
(AT = aacah where acai is the standard deviation of 
the template map). The analysis is done for three differ- 
ent cuts: 20°, 30° and the Tenerife cut (which consists of 
data with 160° < RA < 250°, corresponding to Galactic 
latitudes |6| > 40°). Note that the fits are done jointly 
for M = 2 templates. The DIRBE and Haslam correla- 
tions listed in Table 1 correspond to joint 100/zm— Has 
fits, whereas the R&R numbers correspond to a joint 
lOO/xm-R&R fit. 
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Statistically significant (> 2a) correlations are listed in 
boldface. The two synchrotron templates are found to be 
correlated with the 10 GHz data at low Galactic latitudes 
and this correlation persists with the Has template even at 
higher latitudes, while the 15 GHz data is not correlated 
with these templates even for a 20° cut. The 100/im and 
10 GHz data are correlated only at lower Galactic latitudes 
(20° cut), while the correlation with the 15 GHz data per- 
sists to higher latitudes. The same analysis was carried 
out for the DIRBE 140/xm and 240/mi maps, giving very 
similar results. 

3.2. Latitude Dependence 

To investigate the dependence of the correlation on 
Galactic latitude, we sliced the maps into six regions of 
equal area, each corresponding to a range of latitude 
Figure 2 shows the results for the 100/im map. Note that 
a from the 100/im— 15 GHz correlation is similar in all 
latitude bands, suggesting that the particular ISM prop- 
erties that are responsible for the correlation do not vary 
strongly with latitude. 
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Fig. 2 — Dependence of a on Galactic latitude for the 
100/im— correlated emission at 15 GHz. 



3.3. Spectral Index 

Writing the frequency dependence as 5T oc v" , we 
obtain interesting limits on the spectral index (3. For 
synchrotron radiation, both the Has— 10 GHz and the 
R&R— 10 GHz correlations are consistent with j3 ~ —3 
in all latitude slices. The most favored values between 
408 MHz and 10 GHz values (-3.2 < (3 < -3.4) 
are slightly steeper than the canonical sub-GHz slope of 
-2.7 < < -2.9 (Davies et al. 1998; Platania et al. 
1998). Such a steepening of the spectrum at higher fre- 
quencies is consistent with a steepening of the spectrum 
of cosmic ray electrons at higher energies (Rybicki and 
Lightman 1979). For the emission component correlated 
with the 100/im map, the spectral index between 10 and 
15 GHz is 0.80t;^ for \b\ > 20°, and the spectrum is also 
clearly rising at |6| > 30° (see Table 1). Calibration un- 
certainties are below 5% in the Tenerife data at both 10 
and 15 GHz, and therefore do not affect the conclusions 
that we will draw from this spectral behavior. 



3.4. Effect of Template Correlations 

Equation ^ shows that we can interpret S _1 as mea- 
suring correlations between the various templates. Ta- 
ble 2 shows the dimensionless correlation coefficients r = 
Sr.^S^Sj/)- - 5 , using the 10 GHz pixel weighting, and 
enables a number of conclusions to be drawn: 

1. The three DIRBE maps are, as expected, highly cor- 
related, i.e., they trace the same warm dust emission; 

2. The two synchrotron maps are less well correlated, 
with only r 2 ~ 50% of the variance in the Has map 
being traced by the R&R map at |6| > 20° cut, and 
even less at higher latitudes 1 ; and 

3. All three DIRBE maps are seen to be almost uncor- 
rected with both radio maps, with the common vari- 
ance r 2 being only a few percent. This means that 
S is almost diagonal, and that simply performing 
a separate correlation analysis for each component 
(with N = 1, as was done in e.g., de Oliveira-Costa 
et al. 1997; dOC98) will give almost identical results 
to a joint (N > 1) analysis. We tested this, and in- 
deed obtained results i n go od agreement with those 
presented in subsection 3.1). 



Table 2. — Correlation between foreground templates^ . 
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'"' The upper right triangle of the table uses a |6| > 20° cut 
and the lower left triangle a \b\ > 30° cut. 



3.5. Fake Skies 

We tested our software by analyzing constrained real- 
izations of the CMB and the Tenerife instrument noise. 
Cholesky decomposing the covariance matrix as C = LL T , 
we generated fake skies using y = Lz + Xa, where z is 
a vector of independent Gaussian random variables with 
(z) = and (zz*) = I (the identity matrix), which gives 
(yv T ) - (y)(y>* = L(zz T )L T = LIL T = C. Analyzing 
1000 such realizations, we recovered unbiased estimates a 
with a variance in agreement with equation (0) . 

3.6. Sky Rotations 

Systematic errors or unmodeled foregrounds may add 
non-Gaussian fluctuations to the the data set y, making 
the true error bars larger than suggested by equation (0) . 
To address this issue, we repeated the fit done in equa- 
tion @ after replacing the template map with a set of 
"control" patches selected from different portions in the 
sky. We repeated the analysis with 24 x 2 x 2 = 96 
transformed maps, rotated around the Galactic axis by 
multiples of 15° and/or flipped vertically and/or horizon- 
tally. The correct Haslam template has the highest of all 



1 Systematic effects (such as striping) are present in the synchrotron radio maps and may be partially responsible for the uncorrelated part 
of their signal (Davies et al. 1998). In addition, detector noise in the template maps will of course also reduce their correlation with each other 
and with the Tenerife data, but is in our case negligible compared to the Tenerife noise. 
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96 correlations with the 10 GHz data. Even for a \b\ > 40° 
cut, 95 of the 96 patches (or 99% of them) are less corre- 
lated with the the 10 GHz data than the correct 100/im 
patch. The same test was carried out for other Galactic 
cuts and the 15 GHz data, giving similar results. These re- 
sults show that our correlations are not due to systematic 
errors or chance alignments between the CMB and the var- 
ious template maps. We also applied the rotation test to 
the foreground template correlations in Table 2. Although 
as mentioned above, the dust-synchrotron correlations are 
tiny (r <C 1) and negligible for our purposes, they are still 
statistically significant; for a 40° Galactic cut, the correct 
synchrotron template (Has or R&R) has the highest of all 
96 correlations with the DIRBE map. 

4. CONCLUSIONS 

We have detected high-latitude galactic foreground sig- 
nals in the Tenerife data at high (> 5a) statistical signifi- 
cance, summarized in Figure 3. The synchrotron signal is 
consistent with previous upper limits (Kogut et al. 1996b) 
and detections (dOC98), and provides a useful normaliza- 
tion for modeling microwave foregrounds. We also verify 
that the 15 GHz Tenerife-cut data is useful for CMB cos- 
mology; although the DIRBE-correlated component is sta- 
tistically significant, it is much smaller in amplitude than 
the primary cosmological signal. 

Our results indicate that the DIRBE-correlated signal 
may turn out to be the dominant foreground seen by the 
the MAP satellite. What is its physical origin? A "smok- 
ing gun" indication of spinning dust grains would be a 
turnover in the spectrum, since all Draine and Lazarian 
models show such a feature whereas the antenna temper- 
ature of free-free emission continues to rise with (3 ~ —2.1 
towards lower frequencies. Figure 3 combines our present 
results with previous measurements. The lower panel 
shows that the synchrotron spectrum is well fit by a power 
law v& with j3 ~ —3 (cyan line). In the upper panel, the 
dashed curve shows the spectrum for a best fit linear com- 
bination of vibrational (green) and rotational (yellow) dust 
emission. The vibrational spectrum assumes a 20K dust 
temperature and an emissivity of 2. The rotational spec- 
trum is the warm interstellar medium model of Draine and 
Lazarian (1998) from www.astro.princeton.edu/~draine/. 
A joint three-parameter fit, including free-free emission 
with spectrum i/~ 2 - 15 as a third component (magenta), 
prefers negligible amounts of free-free above 10 GHz. On 
purely physical grounds, free-free emission must of course 
be present at some level. However, if the 10 GHz Tener- 
ife correlation were entirely due to free-free emission, this 
component would explain only about 10% of the vari- 
ance observed at 15 GHz and a percent of the variance 
in the Saskatoon data. The interpretation of such fits 
that directly compare data points from different experi- 



ments is complicated by conversions between different an- 
gular scales (dOC98). However, a direct comparison be- 
tween the 10 and 15 GHz Tenerife points from this work 
is straightforward, since the angular resolution and win- 
dow functions are the same, and gives a spectral index of 
(3 — 0.80^'gg. This is clearly inconsistent with free-free 
emission, and provides evidence that the main culprit is 
indeed spinning dust. 
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Frequency | GHz | 

Fig. 3 — Frequency dependence of DIRBE-correlated emission 
(top) and Haslam-correlated emission (bottom). The DIRBE- 
correlated emission is seen to be approximately fit by a combi- 
nation of spinning and vibrating dust (dashed curve), whereas 
free-free emission alone (falling straight line) cannot explain 
the drop from 15 to 10 GHz. The slope of the Haslam- 
correlated emission is seen to be weekly constrained, fitting 
e.g. a single (3 = —3.26 power law down to 408 MHz (where 
a = 10 6 /iK/K= 1 by definition). Tenerife data is from Table 1 
for a 20° cut. Upper limits are 2-cr. Note that the error bars 
on the 15 GHz point (top) are too small to be visible. 
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